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Abstract. Integral field spectroscopy has been obtained for very few evolved planetary nebulae (PNe). Here we aim 
to study the physical and kinematical characteristics of the unstudied old planetary nebula (PN) PN G342.0-01.7, 
which shows evidence of interaction with its surrounding interstellar medium. We used Integral Field Spectra 
from the Wide Field Spectrograph on the ANU 2.3 m telescope to provide spectroscopy across the whole object 
covering the spectral range 3400-7000 A. We formed narrow-band images to investigate the excitation structure. 
The spectral analysis shows that the object is a distant Peimbert Type I PN of low excitation, formally of 
excitation class of 0.5. The low electron density, high dynamical age, and low surface brightness of the object 
confirm that it is observed fairly late in its evolution. It shows clear evidence for dredge-up of CN-processed 
material characteristic of its class. In addition, the low peculiar velocity of 7 km s _1 shows it to be a member of 
the young disk component of our Galaxy. We further determined an average expansion velocity of V exp = 20.2±1.3 
km s _1 , a local standard of rest radial velocity RVlsr = —27.7 ± 1.7 km s _1 , and a distance of 2.06 ± 0.6 kpc 
for the object. We built a self-consistent photoionisation model for the PNe matching the observed spectrum, 
the H/f luminosity, and the diameter. On the basis of this we derive an effective temperature logT e ff ~ 5.05 
and luminosity 1.85 < logL < 2.25. The temperature is much higher than might have been expected using the 
excitation class, proving that this can be misleading in classifying evolved PNe. PN G342.0-01.7 is in interaction 
with its surrounding interstellar medium through which the object is moving in the south-west direction. This 
interaction drives a slow shock into the outer PN ejecta. A shock model suggests that it only accounts for about 
10% of the total luminosity, but has an important effect on the global spectrum of the PN. 
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1. Introduction 

The planetary nebula PN G342.0-01.7 was discov¬ 
ered and classified as a true planetary nebula in 
the Macquarie/A AO/Strasbourg Ha Planetary Nebula 
Catalogue (MASH I) where it was given the name 
PHR1702-4443 (Parker et al. 2006). Over 1200 true, likely 
and possible new Galactic planetary nebulae were dis¬ 
covered in this AAO/UKST Ha survey of the southern 
Galactic plane. The MASH I & MASH II (Miszalski et al., 
2008) catalogues increase the number of known Galactic 
PNe by ~ 60%, which represents the largest ever incre¬ 
mental sample. According to the morphological scheme 
used in the MASH catalogue, the object was classified as 
an elliptical PN with apparent angular size 90.0"x 63.0". 
Furthermore, it is subclassified as a one-sided enhance¬ 
ment /asymmetry with a well defined “Ear-like” ring struc- 
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ture. The object shows a broken shell with extended dif¬ 
fuse emission seen across a highly filamentary shell struc¬ 
ture (Parker et al. 2006). Weidmann & Diaz (2008) de¬ 
termined the spatial orientation of the object within the 
Galaxy. They measured the position angle PA of the pro¬ 
jected major axis over the sky with respect to the equato¬ 
rial coordinate system EPA = 90°. The EPA was trans¬ 
formed into the galactic system GPA = 142°. This ob¬ 
ject represents a good example of a low surface bright¬ 
ness PN caught at a late stage of its evolution, and pos¬ 
sibly in an early stage of interaction with the surround¬ 
ing ISM. Due to their extreme faintness and large spatial 
extent, such objects have not been spectroscopically stud¬ 
ied in great detail hitherto, although Danehkar, Parker & 
Erolano (2013) have examined one such object, SuWt 2, 
in detail. 

Studying evolved PNe and their central stars provides 
valuable constraints on the evolution from the PN towards 
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the White Dwarf (WD) phase. Whenever these objects 
show signs of interaction with the ISM, they become im¬ 
portant in understanding how the nuclear processed mate¬ 
rials ejected in the PNe are returned to the ISM (Tweedy 
et al. 1995, Tweedy & Kwitter (1996); Xilouris et al. 1996; 
Kerber et al. 1998; Guerrero et al. (1998); Kerber et al. 
2000; Ali et al. 2012). 

The composition of the nuclear processed materials 
ejected into the ISM in the Asymptotic Giant Branch 
(AGB) phase and ionised in the PNe phase are strongly 
dependent on both the initial abundances in the precursor 
star and on the precursor mass, most recently discussed by 
Karakas, 2014 (see also references therein to earlier work). 
In particular, in the more massive, younger population of 
PNe, both Nitrogen and Helium are strongly enhanced by 
the operation of hot-bottom burning during helium shell 
flashes with the subsequent dredge-up of these elements 
into the outer envelope. The PNe with these more massive 
precursors are distinct both kinematically and morpholog¬ 
ically from other PNe, as first recognised by Greig (1971). 
Their morphologies are characterized by the appearance 
of filamentary structure. They also frequently display a 
bi-polar structure. This may have a number of origins. 
These include rotation in the precursor star (e.g. Georgiev 
et al. 2011), binarity of the central star and interaction 
of their stellar winds (Dijkstra, & Speck, 2006; Riera et 
al. 2012), and/or magnetic shaping (Garcia et al.1999). 
Kinematically they belong to the thin-disk population, as 
is to be expected from the relatively short lifetimes of the 
precursor stars. They show a small deviation between the 
velocity determined from the observed local standard of 
rest radial velocity and that determined from the rotation 
curve (Greig 1972, Peimbert and Torres-Peimbert 1983). 
The chemical characteristics of this class of PNe was first 
recognised by Peimbert (1978) and they are now known 
as the Peimbert Type I PNe. 

The evolved Type I PNe, although rather numerous, 
have been very poorly studied. Recently Frew et al (2006) 
discussed two nearby planetary nebulae, similar to the PN 
G342.0-01.7 in the sense they were classified as Type I PNe 
and show low surface brightness, electron density, excita¬ 
tion class, and peculiar velocities. These have the classi¬ 
cal bipolar morphologies and high central temperatures, 
which probably places them in the intermediate age cate¬ 
gory (~ 20000 yr). Manteiga et al. (2004) studied another 
two evolved PNe of low surface brightness towards the 
Galactic bulge. Both show high N/O ratios and He abun¬ 
dances, similar to those observed in Type I disk PNe. A 
large low surface brightness annular PN was detected by 
Pierce et al. (2004), which shows a very low density and 
excitation class with early stage of interaction with the 
ISM. 

The main objective of the present study is to inves¬ 
tigate the global physical and kinematical characteristics 
of the PN G342.0-01.7. For this we have made a mosaic 
of observations covering the whole nebula using the Wide 
Field Spectrograph (WiFeS) - an integral field spectro¬ 
graph mounted on the 2.3m telescope at Siding Spring 


Observatory. The major advantage of such observations 
is that one can extract a global spectrum of the entire 
nebula, avoid contamination by of the nebular emission 
by stars, and derive the global physical and kinematical 
nebular properties such as electron temperature, density, 
chemical abundances, radial and expansion velocities. In 
addition, it is easy to construct an emission line maps for 
the object in chosen emission lines such as [O III], [N II], 
[S II], Ha, and Hf3. 

In this paper, we present our observational data in 
Section [2j discuss the general global properties of PN 
G342.0-01.7 and evidences for its interaction with its lo¬ 
cal ISM in Section [3] and in Section |4] we present the 
results of detailed photoionisation and shock modelling of 
the ejecta. Finally in Section[5]we present our conclusions. 


2. Observations & data reduction 

2.1. ESO observations 

The best overall image of PN G342.0-01.7 has been ac¬ 
quired with the FORS2 instrument (Appenzeller et al. 
1998) mounted on the Cassegrain focus of the Unit 1 
(Antu) of the Very Large Telescope in Chile. PN G342.0- 
01.7 was observed with FORS2 in Period 89 (P.Id.: 089.D- 
0357(A); P.I.: Boffin) as part of series of observations of 
PNe in the MASH catalogue (Parker et al. 2006). 

The extended structure of PN G342.0-01.7 is revealed 
in the Ha+[NH] image (Figure [1]) which is comprised 
of two single 150s exposures acquired on 2012 July 11 
with the “Ha+83” narrow-band filter (Ao = 6563 A, 
A fwhm = 61 A) and downloaded from the ESO archive. 
We reduced both exposures, offset along the North-South 
direction by 1 aremin to cover the gap between the two 
FORS2 CCDs in imaging mode, using the dedicated FORS 
pipeline (v5.0.11) inside the ESO REFLEX environment 
(Freudling et al. 2013). Each four individual CCD images 
(two per exposure) were processed individually and later 
combined into one global, averaged mosaic using a cus¬ 
tom PYTHON script. The data was acquired with 2x2 
binning, so that pixels in the final mosaic have a size of 
0.5 x 0.5 square arcsec, while the seeing at the time of the 
observation was of the order of ~0.8 arcsec. 

As noted above, the effective bandpass of this ESO 
image is A fwhm = 61 A. Thus, all of the [NII] emission 
is included in the bandpass of the filter along with the 
Ha. Since the [NII] emission is dominant in this object, 
the morphology of the PNe in the ESO image is dominated 
by [Nil], 

Figure [T] shows that PN G342.0-01.7 is composed of a 
series of bright filaments oriented along a North-East to 
South-West direction, as well as fainter and shell-like dif¬ 
fuse emission. In particular, we note the presence of diffuse 
emission or an outer envelope - not observed by WiFeS - 
extending towards the North and North-East directions 
(at and around R.A.: 17 ?l 2 rn 5.0 s ; Dec.: -44°42'15.0") with 
respect to the brightest clumps. 
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Fig. 1 . Ha+[NII] emission from PN G342.0-01.7 acquired with F0RS2 on the VLT. The brightest emission clumps 
are oriented along a North-East to South-West axis, with additional diffuse emission forming a shell-like structure. 
Faint whiskers are also extending to the North and North-East from the brightest features. 


2.2. WiFeS observations 

The integral field spectra of the PN G342.0-01.7 were 
obtained over the two nights of March 30-31, 2013 with 
the Wide Field Spectrograph (WiFeS; Dopita et al. 2007, 
2010) mounted on the 2.3-m ANU telescope at Siding 
Spring Observatory. The WiFeS instrument delivers a field 
of view of 25"x 38"at a spatial sampling of either 1.0"x 
0.5"or 1.0"x 1.0", depending on the binning on the CCD. 
For these observations we used the 1.0"x 1.0"option. The 
observations covered the blue spectral range 3400-5700 
A with a low spectral resolution of B ~ 3000 that cor¬ 
responds to a full width at half maximum (FWHM) of 
~ 100 km/s, while in the red spectral range 5700-7000 A 
we used a higher spectral resolution R ~ 7000 correspond¬ 
ing to a FWHM of ~ 42 km/s. The wavelength scale was 
calibrated using the Cu-Ar arc Lamp with 40s exposure at 
the beginning and throughout the night, while flux cali¬ 
bration was performed using the STIS spectrophotometric 
standard stars HD 111980, HD 031128 and HD16O6170- 
These flux standards are claimed to be at the level of ~ 1% 
photometric accuracy. In addition a B-type telluric stan¬ 
dard HIP 38858 was observed to separately correct for the 
OH and H 2 O telluric absorption features in the red. The 

1 Available at : 

www.mso. anu.edu.au/ bessell/FTP/Bohlin2013/GO 12813 .html 


separation of these features by molecular species allows 
for a more accurate telluric correction which accounts for 
night to night variations in the column density of these 
two species. 

All the data cubes were reduced using the PYWIFES H 
data reduction pipeline (Childress et al. 2014) using the 
procedures described therein. In particular, the photomet¬ 
ric calibration of the WiFeS instrument is very reliable, 
since the instrument is based upon image slicing rather 
than optical fibre technology, and the spatial format is 
fixed. The spaxel to spaxel sensitivity variations are cali¬ 
brated using both internal and twilight flat fields, and the 
spectral sensitivity from the multiple standard star ob¬ 
servations. From the sensitivity function derived from the 
three independent standards observed over the two nights 
we estimate the photometric calibration of the instrument 
to be better than 3% in the wavelength range 3700 to 6800 

A. 

PN G342.0-01.7 has an angular size of 90.0" x 63.0" 
square arcsec. Observing the full spatial extent of this ob¬ 
ject required a mosaic of 3 x 3 = 9 WiFeS fields with 1-2 
arcsec spatial overlap between each field. Each individual 
field was observed 2 x 600s = 1200s on-source, simultane¬ 
ously with the red and the blue arm of WiFeS. 

2 http://pywifes.github.io/pipeline/ 
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Fig. 2. The WiFeS emission line map of the PN G342.0-01.7 superimposed on the DSS-2 (red) stellar field. The image 
is a mosaic of the 9 observed positions on the nebula, where each frame has a field of view 25" x 38". The colour 
image is a composite of [Nil] (red) Ha (green) and [OIII] (blue), and reflects the excitation structure of the nebula. 
The horizontal gap results from a pointing offset error between the first and second night’s observations. 


Each frame was reduced individually using pywifes 
into a flux and wavelength calibrated 3-D data cube. The 
final mosaic was assembled using a custom python script 
(available on demand). Because the World Coordinate 
System (WCS) information in the data FITS headers can 
only be trusted to a level of a few arcsec, the Ha and H/3 
integrated flux maps for each frame are constructed and 
used to identify the best spatial shifts between each two 
fields. Given the seeing conditions during the observations 
(1.5"-2"), we round the spatial shifts to the nearest inte¬ 
gers for simplicity, and to avoid superfluous resampling of 
the data. All data cubes were forced on the same wave¬ 
length grid in the data reduction process, so that no shift 
is required along that axis. 

The final alignment of the different fields were checked 
against the Digitized Sky Survey (DSS) 2 red band image 
of the areaH. The final accuracy of the spatial alignment 

3 The DSS-2 image was obtained from the ESO Online 
Digitized Sky Survey: http://archive.eso.org/dss/dss 


of each individual field is of the order of 1 arcsec. In par¬ 
ticular, the DSS-2 image was essential for combining the 
fields acquired during the first and the second night of 
observations, which are separated by a gap of 1-2 arcsec 
as a result of a pointing offset between the two night’s 
observations. 

A composite [Nil] (red) Ha (green) and [OIII] (blue) 
colour image of the mosaic is shown in Figure 2, where 
North is up and East to the left. This image goes much 
deeper and is of higher (spatial) resolution compared to 
that given in the MASH catalogue. The choice of the colors 
is to emphasise the excitation structure of the nebula since 
[O III] is confined to the H + and He + zone, while [N II] 
emission is found only in the H + and He 0 zone of the 
nebula. 

It is apparent that the peak of excitation is displaced 
from the geometrical centre of the nebula, at about RA 
= 17h 02m 05s; Dec = 44° 43’ 31” (J2000). In addition, 
the nebula surface brightness is enhanced around this ex- 
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citation peak. Overall, the displacement of the excitation 
centre towards the SW (determined from the [O III] im¬ 
age) , the enhancement of the surface brightness in the SW 
quadrant, the overall morphology including the faint halo 
extending to the NE seen in the ESO image (Figure [T]) 
are all suggestive of an interaction with the surrounding 
interstellar medium (ISM), with the nebula moving in the 
SW direction relative to its ISM. 

log [N ll]/Ho! 

0.00 0.15 0.30 0.45 0.60 0.75 0.90 1.05 1.20 1.35 1.50 



A R.A. [arcsec] 

Fig. 3. The spaxel to spaxel variation of the [NII] /Ha 
ratio within PN G342.0-01.7. A signal to noise cut of 10 
has been applied. The strongest [NII] emission relative to 
Ha is clearly located in the brightest outer filaments, and 
in the outer parts of the nebular shell. 

The details of the excitation structure are further re¬ 
vealed in a map of the [NII] /Ha ratio, shown in Figure 
3. We see strong [NII] emission relative to Ha located in 
the brightest outer filaments, and in the outer parts of 
the nebular shell, while the lowest values of the [NII] /Ha 
ratio occur towards the centre. The signal to noise in the 
blue part of the spectrum is too low to enable us to pro¬ 
vide similar maps of blue to red line ratios such as Ha/ 
H/3, so we are unable to trace the detailed spaxel to spaxel 
variations in interstellar reddening. 

3. Derived Global Nebular Parameters 

3.1. Nebular spectrum and reddening 

The global spectrum of PN G342.0-01.7 was extracted 
from the mosaiced data cubes of the blue and the red 
spectra. The scaling between the blue and red spectra was 
improved by measuring the total flux in the common spec¬ 
tral range between the blue and red channels (5500-5700 
A), and applying a 2% scale factor to the red spectrum 
to make the two measured fluxes match. The integrated 
nebular fluxes were derived from the mosaiced image af¬ 
ter subtracting the stellar background stars. The resultant 


global spectrum is shown in Figure 4. Note that high in¬ 
terstellar dust reddening, and the residual blue stellar con¬ 
tinuum from faint stars in the field makes it more difficult 
to accurately measure fluxes below A ~ 4200A. 

The line fluxes and their associated uncertainties were 
measured in IRAF splot 0 task. Line fluxes were inte¬ 
grated between two given limits, over a local continuum 
fitted by eye, using multiple Gaussian fitting for the line 
profile. 

In general, the uncertainty of the measured emission 
lines fluxes in the red is smaller than the blue since the 
blue lines are much fainter due to the heavy reddening of 
this object. The first four emission lines in Table 1, have 
measurement uncertainties between 10%-20%. For other 
lines, the uncertainties is less than 6% for lines stronger 
than H/3 and less than 12% for lines weaker than Hp. 
The uncertainties in the form of the reddening function 
increase the measurement errors considerably, particularly 
for the blue lines. We have estimated these errors by using 
various forms of reddening functions from the literature 
with total to selective extinctions in the range 3.0 < R < 
4.6. 

We have used the Nebular Empirical Abundance Tool 
(NEAT Wesson et al. 2012) for plasma diagnostics. 
From the emission lines list, NEAT implemented a stan¬ 
dard techniques to estimate the interstellar reddening, 
electron temperature and density, ionic and elemental 
abundances. The reddening coefficient c(H/3) was deter¬ 
mined from the ratios of hydrogen Balmer lines, in an 
iterative procedure. A preliminary reddening value was 
derived assuming Ha , H(3, Hj line ratios at electron tem¬ 
perature ( T e = 10000 K) and electron density ( N e = 100 
cm" 3 ). NEAT uses a Monte Carlo technique to propa¬ 
gate uncertainties from line flux measurements to elemen¬ 
tal abundances. 

Here, we select the (default) Galactic extinction law of 
Howarth (1983), which is a functional characterization of 
the one from Seaton (1979). We obtained a reddening co¬ 
efficient c(Hp) = 1.81, E(B-V)= 1.24 and a total H/3 flux 
log F(HP) = -12.88. The line ratios [OIII] A5007/A4959 
and [NII] A6584/A6548 both measured to be close to 3.0, 
in good agreement with theoretical predictions of 2.86 and 
3.0 respectively (Storey & Zeippen 2000). Table 1 lists the 
observed and dereddened line strengths of the object. The 
first two columns give the emission lines identification with 
their laboratory wavelengths, while columns 3 and 4 give 
the observed F(A) and the dereddened line strengths 1(A), 
scaled to HP = 100. 


4 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA) under a co¬ 
operative agreement with the National Science Foundation. 

5 The code, documentation, and atomic data are freely avail¬ 
able at http://www.sc.eso.org/~rwesson/codes/neat/ 
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Fig. 4. The blue and red optical integrated spectrum of 
PN G342.0-01.7. Note that residual stellar continuum in 
the field limits the S/N of the line detection below A ~ 
4300. 

Table 1. Observed F( A) and dereddend /(A) line inten¬ 
sities on a scale F(H/3) = 100. The errors on the de- 
reddened fluxes include the uncertainty associated with 
the assumed form of the interstellar reddening function. 


A lab 

Ion 

A (A) 

7(A) 

3726.03 

[O ii] 

34± 6 

100± 25 

3728.82 

[O II] 

59± 8 

172± 35 

4101.74 

H<5 

14.2± 2.7 

29± 7 

4340.47 

Hq 

32.2± 2.0 

54.6± 5.0 

4861.33 

H/3 

100.0± 1.5 

100.0 ± 1.5 

4958.91 

[O in] 

41.0± 1.0 

37.2± 0.8 

5006.84 

[O m] 

129.6± 1.5 

111.7± 2.8 

5200.26 

[N i] 

33.4± 1.4 

23.7± 1.1 

5754.60 

[N n] 

16.8± 1.8 

7.5± 0.9 

5875.66 

He i 

62.5± 2.1 

25.6± 1.5 

6300.30 

[O i] 

177.5± 1.8 

54.9± 5.0 

6363.77 

[O i] 

62.2± 1.2 

18.5± 1.7 

6548.10 

[N ii] 

1108.9± 2.2 

295 ± 25 

6562.77 

Ha 

1120.0± 1.8 

295± 25 

6583.50 

[N ii] 

3427.4± 2.1 

893± 90.0 

6678.16 

He i 

20.6± 1.0 

5.1± 0.6 

6716.44 

[S ii] 

419.4± 1.0 

101.2± 12.5 

6730.82 

[S II] 

298.4± 1.1 

71.4± 9.0 


3.2. Excitation class 

Three methods are known in the literature to determine 
the excitation class (EC) of planetary nebula, which gen¬ 
erally is a good indicator of the temperature of the central 
star (Gurzadyan 1988; Dopita et al. 1992; Reid & Parker 


2010). At low excitation class (up to about 5) the [O III] 
A5007 /H/3 line ratio is the key indicator, while above this 
the strength of the He II emission at 4686 A provides the 
best indication of excitation class of the nebula. This line 
is only present when the central star has T e g > 85000K. 
In G342.0-01.7 the low [O III] A5007 /Hj3 line ratio, the 
high [Nil] /Ha line ratio, and the absence of the Hell 
4686 line would indicate that the object belongs to the 
very low excitation class EC = 0.5, and is also likely to 
show evidence of CN-processed dredged-up material lead¬ 
ing to an enhancement of N and He in the nebula. The 
effective temperature indicated by the excitation class is 
only T e ff ~ 35000K (Dopita et al. 1992). However, as we 
will show below, for evolved PNe such as the one of the 
current study, this naive classification by excitation class 
can be quite misleading. 

3.3. Ionic and elemental abundances 

From the spectral data of Table 1, we can in principle 
derive the electron density from the [S II] AA6716/6731 
line ratio and estimate the temperature from the [N II] 
(A6548 + A6584)/A5754 line ratio. However, the [S II] 
line ratio is found to be close to its low density limit 
N e ([S II]) = 16 cm' 3 , and so does not provide a mean¬ 
ingful constraint on the density. The [N II] lines pro¬ 
vide an estimate of the mean electron temperature of 
T e ([N II]) = 8200 K for the whole PN. 

A first estimate of the ionic and total element abun¬ 
dances of PN G342.0-01.7 was determined by the NEAT. 
Optical recombination lines (ORL) were used to derive 
ionic abundances for helium, while Collisional excitation 
lines (CEL) were used to derive ionic abundances of nitro¬ 
gen, oxygen and sulfur. The total elemental abundances 
were derived applying the ionization correct factors (ICF) 
of Kingsburgh & Barlow (1994), to correct for unseen 
ions. The reddening coefficient, nebular temperature and 
density, ionic and elemental abundances and their asso¬ 
ciated errors are presented in Table 2. The results indi¬ 
cate a pronounced enrichment in both the helium and 
nitrogen abundances (He/H = 0.16 and N/O = 1.0). 
Peimbert (1978) has defined the Type I PNe as the ob¬ 
jects with He/H > 0.14 or N/O > 1. Later on Peimbert & 
Torres-Peimbert (1983) relaxed the definition by consider¬ 
ing those objects with He/H > 0.125 or N/O > 0.5. More 
rigid criteria to define Type I PNe were given by Maciel & 
Quireza (1999), which led to the definition that those ob¬ 
jects satisfy the conditions He/H > 0.125 and N/O > 0.5 
should be classified as Type I PNe. 

The excess in He and N nebular abundances of Type I 
PNe refer to the objects that underwent hot-bottom burn¬ 
ing and dredge-up of CN-processed material into the enve¬ 
lope of the precursor star. This process is restricted to the 
more massive progenitor stars in the intermediate mass 
range, which are therefore associated with the youngest 
Galactic stars found in the galactic thin disk (z < 300 pc, 
Gilmore & Reid 1983). 
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Table 2. Ionic and elemental abundances of PN G342.0- 
01.7. 


Extinction 

c(R3): 

1.81 ±0.02 

Nebular density and temperature 

[S n] density 

16.24 ± 10 

[N n] temperature 

8190 ± 300 

ORL abundances 


He+/H 

0.16 ±4.78 x 10 -3 

He 2+ /H 

- 

He/H 

0.16 ±4.78 x 10 -3 

CEL abundances 


N+/H 

2 55 x iq- 4+ 2 - 79x1 ° “ 

Z.OJ X 1 U -3.62xl0~ 5 

n 2+ /h 

- 

n 3+ /h 

- 

n 4+ /h 

- 

icf(N) 

1 +0-03 

-L.OT_o.03 

N/H 

3 51 x iq- 4+ 4 02x1 ° 5 

0-01 X 1U —3.61 X 10 —5 

O 0 /H 

2 46 x 10 —4+4 23x10 l 
z.no x ru — 3.61 x io —5 

0+/H 
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- 

icf(O) 

1.00 ±0.00 

O/H 

3.49 x 10- 4 ^;“l5 

S+/H 

6.37 x 10- 6 ±?;“-r 

s 2+ /h 

- 

icf(S) 

q /q'T+0.25 
-0.26 

S/H 

„ a r .,„_5+2.79 x10- 6 

2.45 X 10 —2.51 X 10— 6 


The Type I classification of PN G342.0-01.7 on the 
basis of the chemical abundances derived in Table 2 was 
confirmed in light of its kinematic properties. Applying 
the derived PN distance (see section 3.5) the object lies 
at a vertical height z ~ 61 pc below the Galactic Plane. 
In addition it has a small peculiar velocity (|AV|) defined 
as the difference between the observed local standard of 
rest radial velocity and the velocity determined from the 
rotation curve. We find ~ |AV| = 6.9 km s -1 . Maciel & 
Dutra (1992) show that objects belonging to the Galactic 
thin disk should have |AVj < 60 km s -1 . Therefore we 
may conclude that PN G342.0-01.7 belongs to the Galactic 
thin disk consistent with its Type I classification. 

3.4. The search for a central star candidate 

We inspected the digitized UKST (Blue & Red) plates 
from the SuperCOSMOS sky survey, to search for a blue 
central star (CS) candidate. Only one possible CS candi¬ 
date was found which is brighter in the blue plate relative 
to the red plate when compared to the field stars. This 
lies at a = 17:02:4.388, S = -44:43:24.74 (J2000), and is 



Fig. 5. The spectrum of the central star candidate of PN 
G342.0-01.7 extracted from the WiFeS data cube. We con¬ 
clude that this is a field A-Type dwarf star, not related to 
the PNe. 

located very close to the geometric center of the nebula 
(a = 17:02:4.3, S = -44:43:20). Although it lies close to a 
much brighter star, we were able to extract a fairly clean 
spectrum of this object from the WiFeS data cube. This is 
shown in Figure 5. From this spectra we conclude that the 
star is a foreground field A-Type dwarf, and is therefore 
unlikely to be associated with the PN. 

As we will show below, the central star has a very 
high temperature, is very reddened, and therefore would 
remain undetectable in these observations. 

3.5. Distance and H/3 luminosity 

The determination of the distance of a planetary nebula 
is the fundamental parameter which permits an estima¬ 
tion of most nebular characteristics. However, finding a 
reliable distance for the PN is a long-standing problem. 
Generally, two major techniques are used, either individ¬ 
ual estimates based upon, for example, the determination 
of both logT e ff and logc/ of the central star. Failing such 
direct techniques, we must instead rely on statistical meth¬ 
ods. A summary on the applications, assumptions, limita¬ 
tions, and uncertainties of these methods was given in Ali 
et al. (2015) & Frew et al. (2015). 

Many of the statistical methods are based on an em¬ 
pirical relationship between two distance-dependent neb¬ 
ular parameters (e.g. Van de Steene & Zijlstra 1995; 
Stanghellini et al 2008; Ali et al. 2015). All these meth¬ 
ods are rely on the observed 5GHz radio flux. Since there 
is no radio measurement for PN G342.0-01 at 5GHz, we 
have used Equation (6) in Cahn et al. (1992) to convert 
the total H(3 flux into 5GHz radio flux. Using the derived 
5GHz radio flux (30.06 mJy) and the nebular angular ra¬ 
dius (38.25 arcsec, Parker et al. 2006) and following the 
statistical distance scale of Ali et al. (2015), we derived a 
nebular distance of 2100 ± 630 pc. 

Frew et al. (2015) introduced a linear empirical rela¬ 
tionship between the nebular Ha surface brightness and 
nebular physical radii {Sua — r) with a slope between - 
3.3 and -3.8 depending on the calibrated subset used. They 
found that optically thick PNe tend to populate the upper 
boundary of the trend with a mean dispersion of ±28%, 
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while optically thin PNe fall along the lower boundary 
with a mean dispersion of ±18% in the Shoi ~ r plane. 
Applying the definition of optically-thick PNe have rela¬ 
tively strong low excitation features such as [N II], [O II] 
and [S II] lines (Kaler & Jacoby 1989 and Jacoby & Kaler 
1989), and hence the flux ratio of [N II] A6584 /Ha > 1, 
we were followed the optically thick form of Shu — r rela¬ 
tion (Equation 17; Frew et al. 2015). From the spectrum 
analysis, we derived a total Ha flux log F (Ha) = —11.77 
erg cm -2 s . This value was converted into a redden¬ 
ing corrected Ha surface brightness log S (Ha) = —3.55 
erg cm -2 s -1 sr -1 , applying Frew et al. (2015) procedure. 
This surface brightness yields a distance of 2022 ± 566 
pc. We have to take the uncertainty in the calculated dis¬ 
tance with caution, because it is derived from the dis¬ 
persion in the calibration samples used in both distance 
scales. Consequently the uncertainty in the distance could 
be significantly larger. 

Finally we adopt a nebular average distance of 2060 ± 
600 pc. On the basis of this distance, we were derived an 
absolute luminosity of ±h /3 = 4.2 ± 2.5 x 10 33 erg s _1 , and 
a physical dimension of 0.90 x 0.63 pc. 

3.6. Radial and expansion velocities 

We have used the red channel of high resolution spectra 
to estimate the radial velocity. Five emission lines of Ha, 
[N II], and [S II] in the red spectrum were used to find the 
observed radial velocity RV 0 b s of the object. We deter¬ 
mined RV 0 bs = — 57 ± 1.7 km s _1 . This leads to a derived 
heliocentric radial velocity RV h e i = —32.7 ± 1.7 km s _1 . 
The correction to local standard of rest was performed us¬ 
ing the appropriate RV \ ie i of the solar motion. We estimate 
-RVlsr = —27.7 ± 1.7 km s _1 for this object. 

The expansion velocity of a PN is a key parameter in 
determining its evolution. The well defined emission lines 
[N II], [S II] and Ha were used to determine the expansion 
velocity of the PN G342.0-01.7. Following, for example, 
Gieseking et al. (1986): 

Wxp = 0.5 [w 2 obs - < st - 8(In 2 )kT e /m \ 1/2 (1) 

where w obs is the observed FWHM of the measured line, 
Winst is the instrumental FWHM (see Section 2.2), k 
Boltzmann’s constant, T e = 8190K is the electron temper¬ 
ature, and m the atomic mass of measured line. The ob¬ 
served FWHM of each line was determined using the IRAF 
splot task. We determined I4 xp [N II] = 18.8 km s _1 , 
14x P [S II] = 21.2 km s _1 , and F exp [Ha] = 20.6 km s _1 . 
Therefore, we adopt a mean expansion velocity of 14 xp — 
20.2 ± 1.3 km s _1 for this PN. 

3.7. Nebular age 

Provided that we can treat the nebula as being in free 
expansion at V exp = 20.2 km s _1 , the observed size leads 
to a dynamical age of the nebula of ~ 20000 yr. Given 
the evidence presented in Section 13.81 the nebula is most 


Table 3. The derived global parameters for PN G342.0- 
01.7. 


Distance 

2.06 kpc 

Angular size 

90.0"x 63.0" 


Parker et al. (2006) 

Physical size 

0.90 x 0.63 pc 

Volume density 

< 20 cm -3 

Dynamical age 

> 20000 yr 

RVlsr 

-27.7T1.7 km s^ 1 

Vexp 

20.2±1.3 km s _1 

\AV\ 

- 7 km s" 1 

Excitation class 

0.5 

Peimbert class 

Type I 

E(B-V) 

1.24 


likely no longer in free expansion, and therefore these ages 
of the nebular shell should be treated as a lower limit. The 
ionized mass, density and dynamical age of the object all 
support the idea that we are dealing with a considerably 
evolved planetary nebula. A summary of the derived neb¬ 
ular properties is given in Table 3. 

3.8. Evidence for interaction with the ISM 

The PN G342.0-01.7 shows three distinct characteristics of 
an on-going interaction with the surrounding ISM. First, 
the morphology of PN G342.0-01.7 as seen in Figure 1 is 
that of an elliptical ring with a brightness enhancement 
on the southwest side. In addition three straight filaments 
(or stripes) cross the PN in the NE-SW direction, and all 
are brightest in the SW part of the nebula. Such features 
have been seen before in other PNe, e.g. PN G218.9-10.7 
(Ali & Pfleiderer 1999) and NGC 6894 (Soker & Zucker 
1997), where they are also attributed to the interaction 
between the PN and its surrounding ISM. In particular, 
stripes arise due to the stripping of the object’s halo by 
ISM (many examples of such objects were discussed by 
Ali et al. 2012) although Soker & Zucker (1997) speculate 
that such stripes are shaped by the interstellar magnetic 
field due to their orientation. 

Also evident is a faint outer halo, extended to some 
considerable distance in the NE direction. These features 
are consistent with ram-pressure stripping of the PNe shell 
on its trailing edge. 

Second, the excitation centre of the PNe as traced by 
the [O III] A5007 emission in Figure 2 clearly shows a peak 
near RA = 17h 02m 05s; Dec = 44° 43’ 31” (J2000). This 
peak (assumed to be associated with the nearby presence 
of the exciting star) is strongly displaced from the geomet¬ 
rical centre of the main ring of the PNe, which is located 
near RA = 17h 02m 06s; Dec = 44° 43’ 45” (J2000). This 
suggests that the expansion of the nebula has been im¬ 
peded in the SW direction, towards the assumed direction 
of motion of the PNe through its surrounding ISM. The 
brightening of the filaments in the S - SW quadrant of the 
nebula would therefore be ascribed to both an intrinsically 
higher pressure or nebular density, as well as a reduced dis- 
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tance from the exciting star, with consequently a stronger 
radiation field. 

Third, the filaments in the SW quadrant are charac¬ 
terised by an anomalously high [Ol] A6300/Ha ratio. As 
discussed below, we associate this with the presence of a 
slow reverse shock propagating into the PNe ejecta caused 
by the ram-pressure of the motion of the PNe through the 
surrounding ISM. 

4. Theoretical Modelling 

4.1. Photoionisation Model 

We have used the MAPPINGS V code to construct a de¬ 
tailed photoionisation model of PN G342.0-01.7. The 
MAPPINGS V code is a major upgrade of the MAPPINGS IV 
code described in Dopita et al. (2013), and the new 
physics included in this new version will be described in 
Sutherland et al. (2015, in preparation). This code pro¬ 
vides single framework within which to model the emis¬ 
sion line and continuum spectra of equilibrium ionisation 
objects such as HII regions or PNe while at the same time 
being able to do the same for plasmas which are well out 
of collisional or photoionisation equilibrium such as the 
radiative shocks in supernova remnants (SNR) or in the 
Herbig-Haro objects. In the models described here we have 
used the non-LTE model atmospheres described by Rauch 
(2003) taken from the solar abundance ratio grid, and with 
log g = 7 and 8. Where the requested T e g is not computed, 
MAPPINGS V makes a T 4 interpolation between available 
models at each frequency to obtain the ionising spectrum 
of the central star to be used in the modelling. 

Our model aims to match the global spectrum, the size, 
and the absolute H/3 luminosity of the PNe. Since we do 
not know, a priori , the effective temperature of the central 
star, we ran many models of different effective tempera¬ 
ture. Surprisingly, we could not get a satisfactory model 
for any stars with logT e gf < 4.8, although the excitation 
class of EC=0.5 would be expected to yield logT e g ~ 4.5. 
This difference seems to arise as a result of the evolved na¬ 
ture of PNe, and the relatively low density and ionisation 
parameter characterising the photoionised plasma. The 
ionisation parameter can be constrained from the observed 
[O III]/ [O II] ratio, although for these models the varia¬ 
tion with ionisation parameter is relatively weak. We find 
a photon field averaged ionisation parameter log U = —3.4 
in our best-fit model. The central star temperature is con¬ 
strained in the range 4.95 < logT e gf < 5.08. The princi¬ 
ple diagnostics of the effective temperature are the [O III] 
A5007/I7/3 line ratio which increases with T e g, and the ab¬ 
sence of detectable He II A4686 emission, which appears at 
detectable levels in the models when logT e g > 5.1. In our 
best-fit model presented below we adopt logT e ff = 5.05, 
although other models in the temperature range defined 
above are almost equally acceptable. 

The luminosity of the central star is constrained by the 
requirement that the model reproduce the observed H/3 
flux (Lh/ 3 = 4.2 x 10 33 erg s _1 for an assumed distance of 



Fig. 6. The position of the exciting star of PN 
G342.0-01.7 is marked as a red circle as derived 
from photoionisation modelling compared with 
evolutionary tracks of central stars of PNe from 
Vassiliadis & Wood (1993). The gray area repre¬ 
sents the error ellipse defined by the uncertainties 
in the flux and distance estimates, and by the un¬ 
certainties in the theoretical modelling. Note that 
the tracks are marked in terms of the initial stellar 
mass, rather than of the remnant core. 


2.06 kpc. The pressure in the (isobaric) models is adjusted 
until the correct overall dimension of the PNe is achieved. 
This gave log P/k = 5.85cm _3 K. For the best fit model, 
L/Lq = 118 for or an assumed distance of 2.06 kpc. These 
luminosities, the high T e ff of the central star, and the very 
heavy reddening ensure that the central star will remain 
undetectable at visual wavelength. We estimate it to have 
a visual magnitude of ~ 22.0. 

The inferred position of the central star is shown in 
Figure 4 relative to the H-burning evolutionary track of 
Vassiliadis & Wood (1993). On this figure we have also 
plotted isochrones for the evolutionary age dated from 
the time where the central star achieved a temperature 
of logT e ff = 4.0. It is clear that the inferred position of 
the central star of PN G342.0-01.7 is reasonably consis¬ 
tent with the inferred dynamical age of t > 20000 yr. 
However, the uncertainty in the derived stellar parameters 
precludes any estimate of the initial mass of the central 
star. Nonetheless, we can confirm both the evolved nature 
of the object, and its advanced age. 

As described above, the pressure in the (isobaric) 
model is determined by the requirement that it repro¬ 
duce the observed size of the nebula. This is achieved 
for log P/k = 5.85 cm _3 K, for the assumed distance of 
2.06 kpc. From basic Stromgren theory, the nebular pres¬ 
sure required to match the observed diameter scales with 
assumed distance as £> -1 ' 2 . This pressure yields a mean 
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electron density of n e = 26 cm -3 in the ionised plasma. 
The mean electron temperature in the [N II]-enritting zone 
in the model is predicted to be 10000K, appreciably hotter 
than that observed, 8200K, although this latter figure is 
subject to quite large errors, given the faintness of the [N 
II] A5755 line. 

The line flux predictions of our best-fit photoionisa¬ 
tion model are summarised in Table 4 (column 5). This 
confirms the Type I nature of this PNe, with a large over¬ 
abundance of both He and N being derived. In dusty 
PNe models we must be careful to distinguish between 
the gas-phase and total abundances of the elements. The 
gas-phase abundances are what are fitted by a code such 
as NEAT, while the photoionisation model uses the to¬ 
tal abundances of the heavy elements, accounting for the 
fraction of these that are locked up in dust according to a 
standard dust model. The depletion factors of the heavy 
elements onto dust which we use are derived from the for¬ 
mulae of Jenkins, 2009, extended to the other elements on 
the basis of condensation temperatures and/or position on 
the periodic table. To scale the abundance set to include 
dust we must adopt a base Fe depletion in the nebular 
gas, here taken to be logUpe ~ —1.5. 

In our model, we have He/H = 0.155, N/H = 2.3 x 
10 -4 (total), N/H = 2.1 x 10 -4 (gas phase), N/O = 1.42 
(total), N/O = 1.67 (gas phase), O/H = 1.6x 10 -4 (total), 
O/H = 1.3 x 10 -4 (gas phase), and S/H = 7.5 x 10~ 6 
(total and gas phase). These figures should be compared 
with those derived directly from the one-zone spectrum 
analysis in Table 2. Generally speaking, the agreement is 
satisfactory, with the notable exception of O, for which 
the model gives a much lower abundance than the one- 
zone analysis. The cause of this can be found in the [O I] 
zone. The photoionisation model predicts much weaker 
emission than the observations provide. As we will show 
in the next section, this is a symptom of the presence of 
reverse shock emission. 

4.2. Reverse Shock Emission 

The photoionisation model presented above fails to cor¬ 
rectly predict the strength of both the [O I] and [N I] 
emission lines. Both are under-predicted by a factor 4 or 
so. However, in Section [3781 we discussed the evidence that 
the PNe is interacting with its local ISM as a result of its 
motion through this medium. The prime effect of such an 
interaction is to drive a reverse shock through the outer 
layers of the nebular shell. It may also drive a forward 
shock through the ISM, but only if the sound speed in 
this medium is less than the velocity of the interaction. 
Otherwise, there is simply a sub-sonic pressure wave in 
the ISM, which will not be directly observable. Because 
the shell of the PNe is rather dense, the velocity of the 
reverse shock will be low; and of order of the expansion 
velocity. Several 2D and 3D hydrodynamic simulations 
were constructed to investigate how the interaction be¬ 
tween the PNe and ISM affect the outer PN structures, 


see e.g. Villaver et al. (2002, 2003); Muller et al. (2004); 
Wareing et al. (2007). Here we investigate the emission 
properties of such a slow shock, and test whether it is ca¬ 
pable of producing the observed [O I] and [N I] emission. 

Since the PNe is appreciably distorted, we assume that 
the reverse shock velocity is somewhat greater than the 
expansion velocity. For the purpose of the model we take 
V s = 25km s _1 and lognn = 2.0 cm -3 , appropriate to 
the region of cool un-ionised gas ahead of the ionisation 
front in the nebula. The implied ram pressure of the ISM 
is P/k ~ 1.2 x 10 7 cm _3 K. 

A slow shock such as this is characterised by a low 
fractional ionisation of hydrogen at all places through the 
cooling zone. However, as a consequence of the elevated 
post-shock temperature (~ 21000K), we have a large col- 
lisional contribution to the Baimer decrement. In addi¬ 
tion, species such as [O I] and [N I], for which the col¬ 
lision strength for excitation increases with temperature, 
are strongly enhanced. The resultant spectrum predicted 
by the shock model is shown in Table 4 (column 6). Note 
that the more highly ionised species are entirely absent 
in the predicted spectrum, which looks somewhat like an 
extreme low-excitation Herbig-Haro object. 

We find that with only a ~ 13% mixing of shocked 
gas with the photoionised contribution, we can achieve a 
much better fit of the de-reddened observational line in¬ 
tensities (Table 4, column 3) with theory (Tabic 4, column 
4). The shock model provides a surface brightness in H/3 
of logSn /3 = —5.41 erg cm -2 . Since the reverse shock can 
exist only on the leading half of the PNe, we can assume 
that the total area of the shock is 2i tR 2 , where R is the 
nebular radius. On this basis, we estimate the total shock 
luminosity Lh /3 = 2.4 x 10 32 erg s _1 . This represents 6% of 
the total luminosity of PN G342.0-01.7 compared with our 
empirical shock mixing fraction (13%) estimated above. 
Thus we can conclude that the shock model is somewhat 
consistent with the global energetics of the PN. 

5. Conclusions 

In this article, we have discussed the results of the first 
detailed study and optical integral field spectroscopy of 
the PN G342.0-01.7. The spectra cover both the blue and 
red spectral ranges and the data is obtained across a field 
size of 114"x 75". This analysis has provided us with most 
of the physical and kinematical parameters of the object. 
We find that the object is a planetary nebula of appar¬ 
ently low excitation class and in an advanced stage of 
evolution. Plasma diagnostics from both one-zone mod¬ 
els and detailed photoionisation modelling revealed an ex¬ 
cess in both the nitrogen (1.0 < N/O < 1.7) and helium 
(0.15 < He/H < 0.16), placing it in the class of Peimbert 
Type I PNe, with a fairly massive precursor. Kinematical 
measurements show that the PNe is expanding at ~ 20 
km s , and is a member of the thin (young) disk popula¬ 
tion of the Galaxy. A search for the central exciting star 
was unsuccessful. 
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Table 4. Observed line Intensities for PN G342.0-01 compared with our photoionisation and shock models. 




(3) 

(4) 

(5) 

(6) 

A 

Line 

Observed 

PNe ± Shock 

PNe 

Shock 


ID 

(H/3=100) 

(13% shock) 

Model 

Model 

3728 

[O II] 

273 ± 25 

296 

341 

7.9 

4101 

H<5 

29.0 ± 7.0 

23.8 

25.7 

9.9 

4340 

Hy 

54.6 ±5.0 

45.7 

46.7 

39.0 

4861 

H/3 

100 ±1.5 

100.0 

100 

100 

4959 

[O III] 

37.2 ±0.8 

40.6 

47 

0.0 

5007 

[O III] 

112 ±2.8 

117 

136 

0.0 

5199 

[N I] 

23.7 ± 1.1 

23.1 

6.9 

127 

5755 

[N II] 

7.5 ±0.9 

14.4 

16.5 

0.4 

5876 

He I 

25.6 ± 1.5 

22.4 

25.9 

0.2 

6300 

[O I] 

54.9 ±5.0 

54.8 

24.6 

248 

6363 

[O I] 

18.5 ± 1.7 

17.5 

7.9 

79.5 

6548 

[N II] 

295 ± 25 

290.5 

291.5 

15.2 

6563 

Ha 

295 ± 25 

318.5 

291.5 

490 

6584 

[N II] 

893 ± 90 

854.8 

981 

44.8 

6678 

He I 

5.1 ±0.6 

6.5 

7.3 

1.1 

6717 

[S II] 

101.2 ± 12.5 

116 

102.7 

201.5 

6731 

[S II] 

71.4 ±9.0 

81.9 

72.7 

140.2 

Model 

Parameters: 






log L/Lq 



2.07 



logTefl (K) 



5.02 



log g 



7.0 



log P/k (cm -3 K) 



5.85 



Vs (km s -i ) 




25 


log uh (cm -3 ) 




2.0 


We have found clear evidence for an interaction of the 
expanding nebular shell with the surrounding ISM. The 
nebula appears to be moving through this ISM in a south¬ 
east direction. The emission line flux, size, physical den¬ 
sity and global spectrum of PN G342.0-01.7 are matched 
by a model with an effective temperature logT/g- ~ 5.05 
and luminosity 1.85 <; logL <] 2.25. This lies near the 
evolutionary tracks of Vassiliadis & Wood (1993), and is 
consistent with the inferred dynamical age of > 20000 yr. 
The inferred effective temperature of the central star is 
much higher than would be implied by its excitation class 
(EC = 0.5). This is found to be a consequence of the 
evolved state of the nebula and its low ionisation pa¬ 
rameter log U = —3.4. Presumably, this mis-match be¬ 
tween the temperature estimated from the formal excita¬ 
tion class and the true effective temperature of the central 
star would also apply to other evolved PNe of this type. 

A reverse-shock model designed to account for the in¬ 
teraction of the PNe with its ISM suggests that of order 
10% of the flux from PN G342.0-01.7 is shock-excited. 
This fraction is consistent with both the global energet¬ 
ics of the PN, and with the observed nebular spectrum - 
where it is required to explain the anomalous strength of 
[N I] AA5198, 5200 and [O I] AA6300, 6363 emission lines 
relative to H/3. 
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